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Abstract. We present the spectral energy distribution 
(SED) between 10 and 200 /im obtained for the prototype 
Narrow-Line X-Ray Galaxy NGC 7582 with ISOPHOT, 
the photometer on board the Infrared Space Observatory. 
The emission is spatially extended and we separated for 
the first time the nuclear and extranuclear infrared SEDs. 
The nuclear luminosity (Lir ~ 4.5 x 10"'^'' Lq) is dominated 
by cold (T ^ 32 K) dust emission mainly due to star 
formation activity; warm (T ^ 122 K) dust emission is 
also present and is probably related to the active nucleus. 
In addition to a cold component of 30 K, the extranuclear 
SED (LiR - 1.7 X 10^° Lq) shows emission by colder (T - 
17 K) dust: this very cold component comprises 90% of the 
total dust mass of 9.8 x 10^ M©. 
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1. Introduction 

NGC 7582 (z=0.00525) is a highly inclined {i ~ 60°) 
barred spiral galaxy; its blue luminosity corrected for 
galactic and internal extinction is ^ 3.5 x 10^° L© 



(Claussen & Sahai 1992) with a Hubble constant Ho — 
75 km s~^ Mpc"^ It was classified by Mushotzky ( |1982D 
as a Narrow-Line X-ray Galaxy (NLXG) due to its high 
X-ray luminosity, L(2-10 keV) ~ 7x 10® Lq, typical of low- 
luminosity Seyfert 1 galaxies, but with narrow (FWHM < 
200 km s^^) emission lines, more typical of Seyfert 2 or 
starburst galaxies. 

Recent ASCA observations (Schachter et al. 1998, Xue 
et al. 1998) revealed high variability on short {^ 5.5 hours) 
time scales in the hard X-ray range (2-10 keV), typical of 
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Fig. 1. The positions of the circular PHT-P aperture and 
of the ClOO (3x3) and C200 (2x2) arrays are overplot- 
ted on a J-band image of NGC 7582 obtained from the 
SkyView survey analysis system. 



an AGN; in particular, Schachter et al. (1998) report a 
hard X-ray flux drop by 40% in ~ 6 ks, implying a source 
size of '^ 2 X 10^"* cm. In addition, the existence of an 
active nucleus is suggested by the detection of an ioniza- 
tion cone in [O III] A5007 narrow-band images (Storchi- 
Bcrgmann & Bonatto 1991) and by the presence of high 
ionization emission lines (e.g. [Ne V] A3426). Genzel et al. 
(1998) report the detection of both high-ionization (e.g. 
[OIV] X25.9fim, [NV] XU.Sfxm) and low-ionization (e.g. 
[Nell] A12.8/im, [SIII] A18/im) emission lines with the 
ISO short- wavelength grating spectrometer (SWS). The 
low-ionization lines may be produced either by a starburst 



or by an AGN (Radovich et al. 1998). A composite model 



where high-ionization lines are produced in gas ionized 
by the AGN and the low-ionization lines are enhanced in 
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Fig. 2. Intensity distribution in the ClOO array after back- 
ground subtraction. Intensities in the border pixels have 
been normahzed to the range 1-10, with darker colours in- 
dicating increasing intensities; the intensity in the central 
pixel is outside the scale. 
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Fig. 3. Intensity distribution produced by a pointlike 
source in the ClOO array. The intensity in the central pixel 
and the intensity scale are the same as those in the 60/im 
filter for NGC 7582, displayed in Fig. 2. Left: observations 
of the star HR 5340; right: theoretical monochromatic 
footprint as given in Table 2. 



starburst regions, is also possible and would explain the 
[O II I] A5007/H^ ratio of ~ 3 (Storchi-Bergmann et al. 
1995| ), lower than usually obser ved in AGNs. According to 



the 'unified model' (Antonucci 1993 ), broad emission lines 
are extinguished in Seyfert 2 galaxies by a dusty torus lo- 
cated at < 100 pc from the nucleus, but they should be 
scattered into our line of sight by free electrons around the 
torus and should be observed in polarized light. However, 
Heisler et al. ( |1997| ) did not detect broad emission lines 
in the polarized spectrum of NGC 7582: they attributed 
the non-detection of polarized broad lines to the fact that 
our line of sight passes through an amount of dust in the 
torus which is larger than in Seyfert 2s with polarized 
broad lines. 

The presence of strong star formation activity in the 
central kpc was concluded from other indicators. Morris 
et al. (1985) showed that in the inner kpc the Ha emis- 



Table 1. Pointing of the observations. The on-position 
was centered on the object, the off-position was used for 
background subtraction. 



Name 



Date 



NGC 7582 on 
NGC 7582 off 



24 Dec 97 
24 Dec 97 



RA (2000) 
23" 18"" 23^:5 
23'n7'"55n 



DEC (2000) 



-42°22'14" 
-42°20'41" 



sion suggests a rotating disk of HII regions in the plane of 
the galaxy. Recently Aretxaga et al. (1999) reported the 
appearance of broad (FWHM ^ 10000 km s^^) permitted 
lines in the period July - October 1998: they explained 
it as the result of supernova explosions in the circumnu- 
clear starburst, rather than as a change of the reddening 
in the torus which would allow to see the inner nuclear 
regions. The non-thermal radio emission probably origi- 
nates in many supernova remnants associated with this 
starburst, as found by Forbes & Norris ( |1998|) . In the in- 
frared, the 8-13 /xm emission is characterized by the pres- 
ence of polycyclic aromatic hydrocarbon (PAH ) features 



1982| , Roche et 



as reveal ed by ground-based (Frogel et al 
al. 1984) and, more recently, by ISO SWS observations 
(Genzel et al. 1998 ). Their strength is much more typi- 
cal of star bursts rather than of AGNs: in fact, Genzel et 
al. (1998) found that the typical strength of the 7.7/im 
PAH feature, defined in their Table 1, is 0.04 in AGNs 
and 3.6 in starbursts; they find a value of 2.5 in NGC 
7582. Multi-aperture observations performed by Frogel et 
al. ( |1982D showed that the 10 /im emission comes from 
a region smaller than 8", but its origin is still unclear. 
Using the FIR to X-ray luminosity correlation for normal 



and starburst galaxies. Turner et al. (1997) found that the 



maximum contribution of a starburst to the 0.5-4.5 keV 
emission is ~ 34%. According to Xue et al. (199S), the soft 
X-ray emission (0.5 - 2 keV) is the sum of scattered emis- 
sion from the nucleus and emission from the starburst, the 
latter being ^ 10-20%. In their Hubble Space Teles cope 
imaging survey of nearby AGNs, Malkan et al. (1998) de- 
tected dust lanes running across the nucleus of NGC 7582, 
at a distance of some hundreds of parsecs: this dust com- 
ponent external to the torus may contribute to obliterate 
the broad polarized lines whereas the scattered X-ray pho- 
tons may still be observed. 

There is therefore strong evidence that NGC 7582 con- 
tains a 'buried' AGN with the active nucleus residing 
in a dusty environment and coexisting with a circumnu- 
clear starburst. Schachter et al. ( 1998 ) suggested that all 
NLXGs are obscured AGNs. Because of its prototype na- 
ture and brightness we observ ed the infrared SED of this 
object with ISO (Kessler et al. 1996) to determine the rel- 
ative contribution of the various components to the total 
energy budget. 
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Table 2. Footprint values adopted for the spatial de- 
composition of ClOO (Herbstmeier, private communica- 
tion) and C200 (Klaas et al., in preparation) data. For 
ClOO these are calculated values from a convolution of 
the monochromatic point spread function with the ClOO 
pixel field size. For C200 they are empirically derived from 
a comparison of measurements of point sources centered 
on the arrays with raster measurements centered on each 
pixel. 



ClOO Detector 








Pixel 
1 


60/im 
0.0087 


65fj,ra 
0.0132 


80^m 
0.0178 


lOO^m 
0.0228 


2 


0.0429 


0.0423 


0.0379 


0.0385 


3 


0.0087 


0.0123 


0.0176 


0.0228 


4 


0.0440 


0.0434 


0.0377 


0.0382 


5 


0.6625 


0.6587 


0.6416 


0.5561 


6 


0.0440 


0.0434 


0.0377 


0.0382 


7 


0.0087 


0.0123 


0.0176 


0.0228 


8 


0.0429 


0.0423 


0.0379 


0.0385 


9 


0.0087 


0.0132 


0.0178 


0.0228 


C200 Detector 








Pixel 
1-4 


120/im 
0.1435 


150/im 
0.1670 


180/im 
0.1875 


200/im 
0.1875 



2. Observations and Data Reduction 

NGC 7582 was observed with ISOPHOT (Lemke et al. 



1996) in staring mode with the multi-filter AOTs (As- 



tronomical Observation Templ ates) PHT03 (PHT-P) and 
PHT22 (PHT-C) (Klaas et al. |1994[ ) on the target and on 
one off-position (Table 0). The aperture size used with 
PHT-P was 99". On-source integration times were 64s 
with PHT-P and 32s with PHT-C; the same times were 
used for the off-source measurement. The data reduction 
was performed using the ISOPHOT Interactive Analysis 
tool (PIA, version 7.2) together with the calibration data 
set V 4.0 (Laureijs et al. 1998): corrections were made for 
non-linearity effects of the electronics, disturbancies by 
cosmic rays (deglitching) and signal dependence on the 
reset interval time. The flux calibration is based on mea- 
surements with the thermal fine calibration sources (FCS) 
on board. There was one FCS measurement per detector in 
the 12 /zm, 25 /zm, 100 /Ltm and 150 /zm filters . The photo- 
metric accuracies for this observing mode and brightness 
range are ~ 10% in the absolute calibration (FCS in the 
same filter) and < 30% in the r elativ e calibration (no FCS 
in this filter), see Klaas et al. (1998). 



3. Results and Discussion 

3.1. Spatial decomposition 

The redshift of NGC 7582 gives a projected linear scale of 
102 pc/arcsec. The diameter at the 25 m ag arcsec"^ blue 
isophote is 5' (De Vaucouleurs et al. 1991 ), which is larger 
than both the ClOO and C200 detectors (-- 2f2 and 3', 



respectively, see Fig. |l|). Indeed, a first hint that infrared 
emission is extended was found by Sanders et al. (1995) 



in the IRAS Bright Galaxy Sample data reprocessed using 
the HIRES algorithm. They report extended emission in 
NGC 7582 at 12 and 25 ^m (size > 0.'75) and 60 ^m (size 
> 1.'75), but could not resolve it at 100 /zm (size < 2f8). 
The intensity distribution in the ClOO array after back- 
ground subtraction is displayed in Fig. 0. For comparison 
we show in Fig. the patterns expected from a pointlike 
source. The /e/t pattern is that observediov a star, the right 
pattern was generated using the calculated monochromatic 
footprint values of Table g; in both cases the intensity in 
the central pixel was matched to that found in NGC 7582 
(60/um filter). The two distributions are slightly different 
since the observed one includes effects by small telescope 
pointing offsets and chromatic effects within the width of 
the bandpass of the 60/im filter. It can be seen that the 
intensity distribution in NGC 7582 is not that expected 
from a pointlike source. Pixels 2 and 8, which are aligned 
along the major axis of the galaxy, clearly show an excess 
emission; the values in the other border pixels are closer 
to those given by the intensity distribution of a pointlike 
source. 

In order to disentangle the extended and pointlikc 
components, we proceeded as follows. We assumed that 
(i) the emission is composed of a compact central source 
plus extended emission and (ii) the array is centered on 
the compact source. The total intensity observed in each 
pixel is then given by: 



li = .fJn 



L 



(1) 



where /i is the fraction of the intensity from the point- 
like source falling on each pixel (Table y); Iimc is the in- 
tensity emitted by the pointlike nuclear source, Ie,i is the 
extended emission from the extranuclear regions detected 
by each pixel. The total extranuclear emission is defined 
as Je = J2Zi Ie,i (ClOO: np=9, C200: np=4): we shaU caU 
it the disk emission hereafter. For ClOO we made the fur- 
ther assumption that /e,5 is the same as in pixels 2 and 
8, /c,5 = (^,2 + ^e.s) /2; this seems plausible from the ob- 
served intensity distribution on the array (Fig. ||). 

For C200 this approach was not feasible due to the lack 
of a central pixel: we extrapolated the nuclear flux (/iJuc) 
from the blackbody fit of the decomposed nuclear ClOO 
values. The extranuclear component was then obtained 
using Eq. |l| with the footprint scahng for C200 (Table ||): 



T ■ — I- — fl* 



(2) 



An estimate of the uncertainties introduced by this 
procedure was computed by propagating the statistical 
uncertainty of the intensities in each pixel, which is typi- 



cally Ali/Ii < 5% (Laurejis & Klaas |l99g). The compari- 
son (see Table y) with the photometric accuracies quoted 
in Sect, y shows however, that the final uncertainties are 
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Fig. 4. The radio to X-ray spectral energy distribution of NGC 7582. Diamonds: ISOPHOT data (nucleus+disk) ; 
triangles: photometric data fro m the NASA-IPAC extragalactic database (NED), the SX (0.5-2 keV) and HX (2-10 
keV) fluxes are from Xue et al. (1998). Uncertainties are plotted if available. 



dominated by systematic errors in background subtrac- 
tion and relative filter-to-filter calibration which affect the 
total measurement, rather than the individual pixel mea- 
surements. 



3.2. Spectral Energy Distribution 

In Fig. I we show the SED of NGC 7582 from the radio 
to X-ray frequencies, including the new ISOPHOT data 
points. The radio to near-IR SEDs are similar in Seyfert 



and starburst galaxies (Schmitt et al. 1997). The peak of 
the SED is located in the far-infrared and is likely due 
to the dust heated by the starburst. The rise in the blue 
wavelengths of the visual range can be either due to the 



AGN or to the starburst. Schmitt et al. (1999) found that 
nearly 50% of the flux at A5870 A is produced by a popu- 
lation of stars with an age < 100 Myr, and that the con- 
tribution from the nuclear featureless continuum is only 
< 5%; the remaining 45% comes from stars older than 1 
Gyr. The ultraviolet emission was detected by lUE in a 
10" X 20" aperture (linear projected size of 1 kpc x2 kpc). 



Heckman et al. (1995) suggested that the steep ultravio- 
let continuum, fx oc A^'^, may be produced by a very 
dust y sta rburst. This was later confirmed by Bonatto et 
al. ( 1998 ), who reproduced the ultraviolet continuum of 
a sample of nearby spiral galaxies by stellar population 
synthesis and concluded that in the case of NGC 7582 the 
UV light is dominated by recent star formation with an 
age t < 500 Myr in a reddened nuclear starburst. The 
presence of the active nucleus becomes evident only in the 
hard X-ray domain. 

The ISOPHOT spectral energy distribution and its de- 
composition into several spatial and temperature compo- 
nents is presented in Fig. |^. We fitted them with a modi- 
fied blackbody, i.e. a blackbody with a wavelength depen- 



dent emissivity (see e.g. Hildebrand 1983 ); we adopted an 
emissivity ex X~^. Colour correction was performed in an 
iterative way, namely a first value of the temperature was 
adopted and the corresponding colour correction factors 
applied: the resulting SED was fitted with modified black- 
bodies using a non-linear least squares algorithm giving 
new values of the temperature. The whole procedure was 
then repeated until convergence. The fiuxes derived are 
listed in Table 0, the temperature components are given 
in Table |. The fiuxes measured in the PHT-P filters (< 25 
fiTo) with the 99" aperture are in good agreement with the 
values of ground-based observations (Roche et al. 1984) 



with smaller apertures (8"): we therefore conclude that 
an extranuclear component is negligible. The 12 and 25 
fim ISO fluxes are also in good agreement with the peak 
IRAS intensities given by Sanders et al. (1995) with an 
angular resolution (width at 25% peak intensity ~ 94" 
and 70", respectively) comparable to that of the PHT-P 
aperture. 

The results of the decomposition of the nuclear and 
disk SEDs are: 



1. Nucleus - The SED between 10 and 25Aim (PHT-P) is 
well fitted by a warm (T = 122 K) component which is 
typical of Seyfert galaxies (Perez Garcia et al. 1998), 



where dust in the torus is heated by the active nucleus. 
The deviations of the measured data from the fit are 
within the general calibration accuracy, A < 15%, with 
the exception of the 12.8 fiia and 15 /im values (A ~ 
30%). The contribution of the [Nell] A12.8 ^m line 
cannot account for the excess in the 12.8 /im band: in 
fact the intensity given by Genzel et al. ( |1998 ) amounts 
only to ^ 3% of the integrated fiux in this band. At 
A > eO^im (ClOO) we see a wefl defined cold (T - 32 
K) component. Since we have iS',y(25)/S',y(60) = 0.13, 
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Fig. 5. Decomposition of the SED of NGC 7582. Filled symbols: nuclear fluxes; open symbols: disk fluxes. Squares: 
PHT-P; circles: ClOO; triangles: C200 (the 150 ^m disk and nucle ar po ints have been slightly shifted in wavelength 
to avoid overlapping). Diamonds: IRAS data from Sanders et al. ( 1995 ). Asterisks: ground-based observations from 
Roche et al. (1984). The black body fits are indicated as follows. Nucleus: dashed-dotted lines (PHT-P and PHT-C 
components); the thick line gives the sum of the two components. Disk: dashed (ClOO) and dotted (C200) line; the 
thin line shows the sum of the two components. The lower panel displays the deviation of the measured ISOPHOT 
photometric data from the fit. 



whereas in Sey fert galaxies 5'^(25)/S',.(60) > 0.2 (Mi- 
ley et al. 1985| ), we conclude that at these wavelengths 
the emission is dominated by dust heated by the star- 
burst. According to the colour-colour diagrams by Do- 
pita et al. ( |199S| ), the FIR colours of NGC 7582 may be 
interpreted as the composition of substantial circum- 
nuclear star formation emission plus emission from an 
active nucleus observed through a dusty environment 
(''"12/im > 3): the dust may be located both in the ac- 
cretion torus and in the circumnuclear starburst. 
2. Disk - The signal in ClOO is dominated by a com- 
ponent whose temperature (T ^ 30 K) is similar to 
that found in the nucleus. The signal in C200 reveals 
a colder component, T ~ 17 K; due to absence of data 
points longward of 200 /Ltm, the fit is not very well 
constrained and has an uncertainty AT ^^ 3 K. This 
temperature is typical of dust heated by the interstel- 
lar radiation field; the presence of very cold dust which 
dominates the emission in the outer regions of the disk 
has also been found in the analysis of ISOPHOT data 
both for active (Perez Garcia et al. 1998) and normal 
(Kriigel et al. |1998| , Alton et al. |1998| ) galaxies. 



3.3. Luminosities and dust masses 

Luminosities between 1 and 1000 /im derived for the dif- 
ferent components are compiled in Table H. The ratio of 
the nuclear starburst t o disk luminosity is 2^^. Rowan- 
Robinson & Crawford ( 198E ) obtained a value of 2.1 fit- 
ting the IRAS values with starburst and disk emission 
models: this provides a reliability check for the decompo- 
sition based on the spatial resolution we obtained. 

The spatial and spectral decomposition allows an es- 
timate which is more reliable than with IRAS of the 
star formation rates and dust masses (Table 0) . Star for- 
mation rates were derived for the cold (30 K) nuclear 
and disk components only, according to Kennicutt ( 1998 ); 
dust m asses were computed according to Klaas & Elsasser 
(1993), assuming dust grain properties as in Hildebrand 
(1983) for the emissivity adopted here (ex A^^): 



SFR = 1.7 X 10-i°LiR,/io [Mo yr^i], 
Md = 7.9 X 10-5(rK/40)-6 Lin/LQ [Mq]. 



(3) 

(4) 



The nuclear star formation rate of ^^ 6M yr ^ is typ- 



ical of bright IR galaxies (Kennicutt 199S ) , and confirms 
the presence of a moderate star formation activity in the 



6 M. Radovich et al: The 10-200 ^m SED of NGC 7582 

Table 3. Central wavelengths and spectral resolutions of the ISOPHOT filter bands, total fluxes (not colour corrected) 
and fluxes obtained after separation of the nuclear and disk components, blackbody fitting and colour correction; 
statistical uncertainties for the decomposed components (in italics) and photometric accuracies are also given. We 
show for comparison the IRAS fluxes from Sanders et al. ( 1995| ); a proper colour correction is applied only to the peak 
and extended values. Fluxes and uncertainties are in Jy. 













ISOPHOT 










IRAS 






Xc 


A. 
AA„ 


Total 


phot 


Nucleus 


Stat 


phot 


Disk 


Stat 


phot 


Total 




Peak 


Extended 


PI Detector 


























10 


5.6 


0.6 


±0.2 


0.6 




±0.2 
















12 


5.6 


1.5 


±0.2 


1.8 




±0.2 








2.3 


±0.1 


2.1 


0.9 


12.8 


5.3 


2.8 


±0.8 


2.8 




±0.8 
















15 


1.8 


2.6 


±0.8 


2.7 




±0.8 
















P2 Detector 


























20 


2.3 


6 


±2 


6 




±2 
















25 


2.6 


7.2 


±0.7 


5.8 




±0.6 








7.5 


±0.4 


5.1 


0.7 


ClOO Detector 


























60 


2.5 


57 


±17 


42 


± 2 


±13 


11 


±1 


±3 


52 


±3 


54 


6 


65 


1.1 


59 


±18 


53 


± 3 


±16 


16 


±1 


±5 










80 


1.7 


81 


±24 


64 


±4 


±19 


28 


±2 


±9 










100 


2.2 


85 


±9 


62 


±4 


±6 


29 


±2 


±3 


78 


±5 






C200 Detector 


























120 


2.4 


90 


±27 


52 


±3 


±16 


39 


±4 


±12 










150 


1.9 


75 


±8 


36 


±2 


±4 


39 


±4 


±4 










180 


2.6 


59 


±18 


25 


±1 


±7 


32 


±3 


±10 










200 


3.6 


50 


±15 


19 


±1 


±6 


30 


±2 


±9 











Table 4. Aperture diameters and derived projected linear sizes (radii), infrared luminosities, dust masses and star 
formation rates derived for the different components; uncertainties given by the blackbody fit are displayed. 



Aperture 


Radius 


T 


L(l-1000 lira) 


Mdust 


SFR 


n 


[kpc] 


[K] 


[i0] 


[Mq\ 


[^0 yi-"'] 


Nucleus 












99 


< 5 


122±5 


(12lD X 10' 


(111?) X 10^ 


- 


43 


< 2 


32±0.3 


(3311) X 10« 


(93I5) X 10^ 


5.81":^ 


Disk 












138 


< 7 


30±2 


(11+^) X 10^ 


(5+1) X 10'' 


2t\ 


182 


< 9 


17±3 


i&tf) X 10^ 


(Sta) X 10" 





inner kpc. A lower star formation rate, 
derived for the disk. 



2M, 



o yr 



According to Claussen & Sahai ( I992| ) the atomic and 
molecular hydrogen masses in NGC 7582 are M(HI) = 
4.3 X lO'^ M0 and M(H2) = 4.4 x 10^ M©, respectively. The 
dust to gas mass ratio would be ~ 1/600, if we considered 
only the warm and cold components (nucleus ± disk) de- 
tected by ClOO; this would give a dust mass Md ~ 1.4 x 10^ 
Mq. The addition of the 17 K component revealed by C200 
increases the dust mass by a factor of 7; the dust mass is 
now Md ^ 9.8 X 10^ M©, the dust to gas mass ratio is 
1/90, close to the canonical value. It should be noted that 
due to the strong dependence of the dust mass on the 
temperature, the AT = 3 K uncertainty may change the 
dust mass of the very cold component by a factor of 3. 
A more detailed error calculation shows however that the 



final uncertainty on the dust mass is much lower (see Ta- 
ble 0). The dust mass of the nuclear warm component, 
probably associated with the torus, is Md '^ 3 x lO'^ M©. 

4. Summary 

We obtained ISOPHOT data for the prototype NLXG 
NGC 7582 from 10 to 200 /im and found that its emission 
is extended in the 60 to 200 /im range along the major axis 
of the galaxy. We performed a spatial decomposition into 
a nuclear and a disk component. The SED in the nucleus 
(r < 2 kpc) amounts to 70% of the total IR luminosity 
of 6 X 10^° Lq; it was decomposed into a warm (122 K) 
component, probably related to dust heated by the ac- 
tive nucleus, and a cold (32 K) component coming from 
a circumnuclear starburst. The IR colours are consistent 
with the values expected for an obscured active nucleus 
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plus circumnuclear star formation. We conclude that the 
most likely interpretation of the NLXG nature of NGC 
7582 is that the high X-ray luminosity is produced by the 
active nucleus, which is obscured at longer wavelengths 
by dust located in the circumnuclear regions; the heat- 
ing of the dust by the starburst produces the far-infrared 
emission. The availability of the data longward of 100 /im 
allowed the detection of a very cold (17 K) component in 
the SED of the disk. This very cold component is due to 
dust heated by the interstellar radiation field and it dom- 
inates the dust mass of the galaxy; its inclusion increases 
the dust mass detected now by a factor of 7. 
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